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Abstract: We carry out preliminary numerical study of Sugino's lattice formulation [Q, ||| 
of the two-dimensional M = (2, 2) super Yang-Mills theory (2d M = (2, 2) SYM) with the 
gauge group SU(2). The effect of dynamical fermions is included by re- weighting a quenched 
ensemble by the pfaffian factor. It appears that the complex phase of the pfaffian due to 
lattice artifacts and flat directions of the classical potential are not problematic in Monte 
Carlo simulation. Various one-point supersymmetric Ward-Takahashi (WT) identities are 
examined for lattice spacings up to a = 0.5/ g with the fixed physical lattice size L = 4.0/g, 
where g denotes the gauge coupling constant in two dimensions. WT identities implied 
by an exact fermionic symmetry of the formulation are confirmed in fair accuracy and, 
for most of these identities, the quantum effect of dynamical fermions is clearly observed. 
For WT identities expected only in the continuum limit, the results seem to be consistent 
with the behavior expected from supersymmetry, although we do not see clear distintion 
from the quenched simulation. We measure also the expectation values of renormalized 
gauge-invariant bi-linear operators of scalar fields. 
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1. Introduction 

It will be very exciting if non-perturbative question in supersymmetric gauge theories (such 
as possibility of spontaneous breaking of supersymmetry) can be studied numerically at 
one's will. Despite the great efforts being made towards numerical study of the four- 
dimensional Af = 1 super Yang- Mills theory (4d J\f = 1 SYM) §§— so far no conclusive 
evidence of a restoration of supersymmetry in the continuum limit has been observed. 



For recent reviews on lattice formulation of supersymmetric theories, see refs. [10, 11, 12]. 
Under this situation, to test various ideas, it seems useful to examine lower dimensional 
supersymmetric gauge theories in great detail, which have much simpler ultraviolet (UV) 
structure and for which it is relatively easy to accumulate high statistics in Monte Carlo 
simulation. 

In this paper, we report the results of our small-scale Monte Carlo study of lattice 
formulation of the 2d M = (2,2) SYM, proposed by Sugino M, M. For this and simi- 
lar lower-dimensional supersymmetric gauge theories, many other proposals and studies 



on possible lattice formulation exist |L3J-|25|]. (See also ref. |26[ for studies based on the 



supersymmetric discrete light-cone quantization.) The advantage of the formulation of 
refs. (l|, |2| is that a fermionic symmetry, associated with one of four supercharges of the 
target theory, is manifestly preserved even with finite lattice spacings and finite volume. 
Full supersymmetry is expected to be restored in the continuum limit. Possible disadvan- 
tage of the formulation, on the other hand, is that the pfaffian resulting from the integration 
over fermionic fields is generally complex, 1 although the complex phase is expected to be 



1 To avoid this point is one of motivations of the proposal of ref. ]is[]. 
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irrelevant in the continuum limit, as the corresponding pfaffian in the target theory is real 
and positive semi-definite. 



In our simulation, we include the effect of dynamical fermions by re-weighting. That 
is, in taking a statistical average, a quenched ensemble is re-weighted by the factor of 
pfaffian. With parameters and statistics of our Monte Carlo simulation, it appears that 
the complex phase of the pfaffian and flat directions of the classical potential (which might 
imply subtlety in the integration over scalar fields) are not problematic. The parameters 
of our simulation correspond to lattice spacings up to a = 0.5/ g with the fixed physical 
lattice size L = 4.0/ g, where g denotes the gauge coupling constant in two dimensions. 

In this paper, we mainly study one-point super symmetric bare WT identities. These 



are precisely WT identities numerically analysed by Catterall [22] on the basis of his lattice 
formulation of the 2d J\f = (2,2) SYM In our numerical simulation, WT identities 

implied by the exact fermionic symmetry of the formulation are reproduced in fair accu- 
racy and, for most of these identities, we clearly observe the quantum effect of dynamical 
fermions. For WT identities expected only in the continuum limit, the results seem to 
be consistent with the behavior expected by supersymmetry, although we do not see clear 
distinction from the quenched (i.e., non supersymmetric) simulation. We measure also the 
expectation values of renormalized gauge-invariant bi-linear operators of scalar fields to 
illustrate how this kind of numerical study would be useful. 

In section 2, we briefly review Sugino's formulation of the 2d N = (2,2) SYM mainly to 
fix our notation. Some remarks are made on the continuum limit. In section 3, the results 
of our Monte Carlo simulation are reported. In section 3.1, we explain our simulation 
algorithm and related matters. In section 3.2, one-point WT identities are studied. In 
section 3.3, expectation values of gauge-invariant bi-linear operators of scalar fields are 
studied. Section 4 is devoted to conclusion. Throughout this paper, the gauge group is 
assumed to be SU(iV c ) and our simulation has been done only for SU(2). 



2. Sugino's lattice formulation of the 2d M = (2, 2) SYM 



2.1 Topological field theoretical form of the continuum target theory 



This lattice formulation starts with the fact that the (euclidean) action of the 2d J\f = (2, 2) 



SYM can be written in the form of the topological field theory [27]' 



^continuum = ^ I d 2 X tvl -[(ft, 0] 2 + H 2 - iH<S> + D^D't, 



- \vi<t>, v] ~ x[4>, x] + + ixQ$ + i^D^rX , (2.2) 

where all fields are SU(iV c ) Lie algebra valued and scalar fields (ft and (ft are combinations 
of two real scalar fields, (ft = Xi + iX% and (ft = X% — 1X3, respectively. f> = 2F01 is the 
field strength in two dimensions Fqi = 8qAi — d\Ao + i[Aq, A±\. The covariant derivatives 
Dfj, are defined with respect to the adjoint representation D^cp = d^tp + £[-A^,<p] for any 
field (p. The index \x runs over and 1. Note that, in the above convention, the bosonic 
fields A^, (ft and (ft have the mass dimension 1 and the fermionic fields ipu, x an d V have 
the mass dimension 3/2, because the gauge coupling constant in two dimensions g has the 
mass dimension 1. 

In eq. (|2.2j), Q is a BRST-like transformation in the topological field theory (that is a 
particular linear combination of super-transformations in the original SYM theory) and is 
defined by 

QAft = iftfi, Qiftft = iD^, 

Q(ft = 0, 

QX = H, QH=[<ft, X ], 

Q$ = V, Qr}=[<t>A\- (2-3) 

The salient feature of this transformation is that its square Q 2 is an infinitesimal gauge 
transformation with the transformation parameter (ft. Therefore, Q is nilpotent Q 2 = 
when acting on gauge invariant quantities. Moreover, the action can be expressed as a 
Q-exact form: 

^continuum = J ^ X tr { \v[<P, 4>] ~ «X* + XH ~ m/^-D^ j . (2.4) 

In this form, the Q-invariance of the action is manifest. Then the idea 3 is to construct 
a lattice analogue of the Q transformation such that the nilpotency (up to the lattice 
gauge transformation) holds. Then adopting a lattice action of the structure of eq. ([2.4j), 
Q-invariance can be preserved exactly in lattice theory. 

2 The conventional form of the action of the 2d M — (2,2) SYM, for example, eq. (2.7) of ref. |l9|, is 
reproduced by the following substitution 

a a a 

4> -» -ig^-itf + - ^ +^)r a /2, - -ig^2W - ^ + Oft - i^ a 2 )T a /2, 

a a 

X - - ^ - #i - #2)^72, V - "ij^W + ^ +Ti +i>l)T a /2, (2.1) 

a a 

where T a axe anti-hermitian generators of SU(iV c ) normalized as tv{T a T b } = — (l/2)d a b and the index a 
runs from 1 to N% — 1. 
3 See ref. 0. 
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2.2 Lattice formulation 

We consider two-dimensional square lattice of the one-dimensional physical extent L, 

A = {x G al? | < x M < L) , (2.5) 

where a denotes the lattice spacing. We define also the one-dimensional extent in a lattice 
unit N = L/a. All fields except the gauge potentials are put on sites and, as is conventional 
in lattice gauge theory, the gauge field is expressed by the compact link variables U(x,fj,). 
Periodic boundary conditions on A are assumed on all fields. 

As a lattice counterpart of the fermionic transformation (2.3), we define {jx implies a 
unit vector in the /^-direction) 

QU(x,n) = iip,j,(x)U(x,fj,), 

Qipn(x) = i^^{x)tp^{x) - i [<j>{x) - U (x, n)(f>{x + afi)U(x, , 
Q(j)(x) = 0, 

Q X (x) = H(x), QH(x) = [</>(x),x(x)], 

Q4>{x) = nix), Qr,{x) = [<t>{x)A{x)). (2.6) 

It can be confirmed that Q 2 is in fact an infinitesimal lattice gauge transformation with 
the parameter <f>(x). Thus the nilpotency Q 2 = holds on gauge invariant quantities. The 
lattice action is then defined by an expression analogous to eq. (E 



S = Q* 2 E + ° 2 ^ + + ^2 tr ix(x)H(x)}] , (2.7) 

x&A ^ a 9 ' 



where 



Oi{x) = tr|^(x)[0(x),0(x)]| , (2.8) 

2 (x) = tr , (2.9) 

°^ x ) = ^2 tT ^Y,M X ) $( x ) - U(x,^(x + afi)U{x^)- 1 ) J . (2.10) 

In the above expression, <fe(x) is a lattice analogue of the field strength and is defined from 
the plaquette variables 

U(x, 0, 1) = U(x, Q)U(x + aQ, l)U(x + al, 0) _1 ?7(x, l)" 1 (2.11) 

by 



= : — 1 \! 2 (2-12) 



$(x) 

i-^||i^7(x,o,i 

with 

<$>(x) = -i [U(x, 0, 1) - U(x, 0, l)- 1 ] . (2.13) 
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Finally, the matrix norm in the above expression is defined by 

1/2 



A\ 



tr 



(2.14) 



and the constant e is chosen as (for N c = 2) 

< e < 2V2. (2.15) 

(The meaning of the denominator of eq. ( 2,12[ ) will be explained shortly.) From the Q- 
exact form ( |2.7| ) and the nilpotency of Q, the lattice action is manifestly invariant under 
the Q-transformation fl2,6p . 4 

After the operation of Q, the lattice action becomes 

) , (2-16) 



s = a2 E ( E £ b*(x) + E c ^ x ) + i tr { 

xeA \i=i i=i ° 9 1 



1 



F(x) - -i*TL(a?) 



where we have noted that only the traceless part of &(x), 

$ TI » = - — tr{$(x)} 1, 



(2.17) 



appears in the action, because the auxiliary field H(x) is traceless Each term of the 
action density is given by 



Cbi(x) 



1 



M 7 &tl(x) 



if 1 

a 9 l«=0 



(2.18) 
(2.19) 



x (0(z) - C/ (x, n)4>(x + afi)U(x, n) 1 ) 



(2.20) 



and 



C F i(x) 

C F3 (x) 
C Fi (x) 

Cf6(x) 



atg* 



a 4 g 2 



tr |--r/(x)[(/>(j;),r ? (x)]| , 
tr{-x(x)[0(x),x(a;)]} , 



ti {-ipo(x)ip {x) (4>{x) + U(x,0)^(x + a6)U(x,0)- 1 )} , 
tr {-V>i(x)V>i(x) (^(a;) + (x, l)0(x + al)U(x, l)" 1 ) } , 
tr {ix(x)QHx)} , 

7-2 m -*E fafa) ~ ^ ( x ' + a P-) u ( x > ^r 1 ) \ 

9 { J 



(2.21) 
(2.22) 
(2.23) 
(2.24) 
(2.25) 

(2.26) 



4 Another important property of the present lattice formulation is a manifestly preserved global U(1)r 
symmetry || ||. 



- 5 - 



It is important to keep in mind that all lattice fields in the above expressions are dimen- 
sionless. For comparison of correlation functions with the continuum theory ( |2,2j ), we have 
to rescale all lattice fields (and Q) by appropriate factors of 1/a according to their mass 
dimension. 



With the lattice action (2.16), the expectation value of an operator O is defined as 
usual 

'' d^Oe~ s 



«0» = J — , (2.27) 

d/j.e s 



where the integration measure is defined by (writing (f>(x) = X 2 (x) + iX^(x) and <fi(x) = 
X 2 (x) - iX 3 (x)) 

d ^= II ( Il dU ( x ^) ] fl aX2(x)dX^x)dH a (x) I n<¥» ] d X a (x)dr ] a {x) 

(2.28) 

N 2 — 1 

in terms of color components of fields, <p(x) = — i J^a=i ^{x)^, where T a are anti- 
hermitian generators of SU(-/V C ) (normalized as tr{T a T b } = —(l/2)5 a b). dU(x,ii) is the 
standard Haar measure. Note that the integration over the auxiliary field H{x) is gaussian 
and can be done readily. The invariance of this measure under the Q-transformation is 
noted in the last reference of ref. flBfl . 

The denominator of eq. ( 2.12j ) needs an explanation. Without that factor, the lat- 



tice action for the gauge field is the "double-winding plaquette type" [g9| and the action 
possesses many degenerate minima which have no continuum counterpart. Due to the de- 
nominator of eq. ( 2,12j ), the action (2.16) diverges as ||1 — U(x, 0, 1)|| ->eata certain site x. 



Precisely speaking, the above construction of the action is applied only for configurations 
with 

||1 - J7(a?, 0,1)|| < e, forVxeA, (2.29) 
and, otherwise, i.e., if there exists x £ A such that ||1 — U(x,0, 1)|| > e, we set 

S = +oo. (2.30) 

In this way, the domain of functional integral ( |2.27| ) is effectively restricted to the space 
specified by the condition ( p.29| ) 5 and, by setting the parameter e in the range fl2.15| ), it can 
be shown that the unique physical minimum of the action (up to gauge transformations) 
is singled out. This procedure to solve the problem of degenerate minima does not break 
the Q-symmetry. See ref. Q for careful discussion on these points. 

With the above construction, one fermionic symmetry Q is manifestly preserved on the 
lattice. The price to pay is that the pfaffian, resulting from the integration over fermionic 
fields, is generally complex 6 and this could be disadvantage in Monte Carlo simulation. We 
will see below that, however, this point appears to be not problematic, at least with the 
parameters in our numerical study. 

5 This is the so-called admissibility condition considered in a different context ptj. 

6 In the target continuum theory, the corresponding pfaffian is real and positive semi-definite. 
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2.3 Continuum limit 

In the present two-dimensional super-renormalizable system, all dimension-ful quantities 
can be measured by taking the gauge coupling constant g, which has the mass dimension 1, 
as a unit (in this sense, g is analogous to the A-parameter in QCD). 7 The continuum limit is 
defined by the limit a — > 0, while g and L, the physical extent of the two-dimensional space, 
are kept fixed. In refs. |l5|, U, the restoration of full set of supersymmetry in this continuum 
limit was argued on the basis of the loop expansion and power counting. 8 More precisely, 
this argument shows that the 1PI effective action for elementary fields is supersymmetric 
in the continuum limit. Note that the argument of refs. [15| , [l| says nothing about possible 
supersymmetry breaking in correlation functions that contain composite fields. 

Now, in numerical study, it is convenient to define the dimensionless gauge coupling 
constant by 

P 1 



2N C a 2 g 2 ' 



(2.31) 



that is simply the over-all common coefficient of the lattice action (2.16). Clearly, (3 goes 
infinity in the continuum limit. In terms of (3, the lattice spacing in a unit of the gauge 
coupling constant g is given by 

a = ff? (2 ' 32) 

and, correspondingly, the one-dimensional physical extent of the lattice is 




L = aN = a —-N-, (2.33) 



where iV is the one-dimensional size in a lattice unit. 

As already noted, all fields on the lattice must be rescaled by appropriate factors of 1/a, 
for comparison with the continuum theory ( |2.2| ) . All bosonic fields in the continuum theory 
(except the auxiliary field), which have the mass dimension 1, are related to the lattice 
fields by 

i nr~ 

^continuum (%) = ~<f(x) = \ -^T9<P{x) ( 2 -34) 
a V 2iv c 

and the correlation functions are measured in a unit of g. Similarly, fermionic fields are 
related as 

^continuum fa) = ^TJ^fa) = [qIT J # V>fa)- (2.35) 

Note that, in the continuum limit, these rescalings amplify the correlation functions on the 
lattice. 



7 Recall also that, in the present system, there is no non-trivial coupling constant renormalization nor 
wave function renormalization. Only mass terms of bosonic fields may be renormalized (ignoring gauge 
symmetry and supersymmetry). 

8 Strictly speaking, this argument as it stands holds for the limit a — » with the fixed number of lattice 
points N = L/a (thus the physical lattice size goes to zero L = aN — > 0). The argument, however, can 
slightly be modified to show a restoration of supersymmetry in the present (L fixed) continuum limit, to 
all orders of the loop expansion. 
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3. Monte Carlo study 



3.1 Algorithm, simulation code and statistics 

In supersymmetric theories, the quantum effect of fermions is vital and the quenched ap- 
proximation is almost meaningless. Even in the present two-dimensional system, a treat- 
ment of dynamical fermions can be non-trivial and costly. The 2d N = (2, 2) SYM can 
be obtained by dimensional reduction of the 4d N = 1 SYM in which the fermion field is 
a Majorana spinor instead of Dirac. Thus the pfaffian of the Dirac operator, instead of 
the determinant, naturally appears. In a sense, we have to treat an Nf = 1/2 system. To 
compute the pseudo-fermion force in the hybrid Monte Carlo algorithm, one then has to 
implement the fourth-root of D^D, where D is a lattice Dirac operator. Moreover, this 
fermion must be massless (at least in the continuum limit). Thus the numerical simulation 
of the 4d SYM is quite demanding. 

In two dimensions, on the other hand, it should be relatively easy to accumulate high 
statistics compared to four dimensions. Taking these things into consideration, here we 
adopt a (somewhat brute force) re-weighting method. 9 That is, we prepare configurations 
with the statistical weight e~ SB , where 5b is the lattice action ( |2.16| ) with all fermion fields 
are removed. This is a quenched ensemble. Writing the expectation value in this purely 
bosonic system by 10 

/ d/i B Oe- 5B 

(O) = J — , (3.1) 

d/^B e~ SB 



the true expectation value is evaluated by re-weighting configurations by the factor of 
pfaffian 11 



where D is the lattice Dirac operator appeared in the action ( |2.16|) , Mathematically, 
this definition is equivalent to the original one ( 2.27] ). Practically, however, we have only 



a limited number of configurations and there may exist the overlap problem. That is, 
distribution of configurations favored by the the quenched weight e~ SB may not have a 
sufficient overlap with that of configurations really important in the original system. So 
we need many configurations to reproduce the true expectation values in the original un- 
quenched system. 

We developed a C++ code of the hybrid Monte Carlo simulation with the action Sb 
by using a library due to Massimo Di Pierro, the FermiQCD/MDP |33]| . For each con- 
figuration, we compute the inverse (i.e., the fermion propagator) and the determinant of 

9 Hidenori Fukaya suggested this method to me. For application in two dimensions, see, for example, 
ref. §jj. 

10 When the operator O contains fermionic fields, they are contracted by fermion propagators in the 
presence of bosonic fields. 

11 Note that, in this method, ((C)) is evaluated by a ratio of two averages over an ensemble. This means 
that ((C)) is not the primary quantity ]32j and care is needed to estimate the statistical error in ({0}). We 
used the jackknife analysis to estimate the average and the statistical error for ((C)). I would like to thank 
Issaku Kanamori for clarifying discussion on this point. 
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N 


8 


7 


6 


5 


4 


P 


16.0 


12.25 


9.0 


6.25 


4.0 


number of configs. 


1,000 


10,000 


10,000 


10,000 


10,000 


ag 


0.5 


0.571428 


0.666666 


0.8 


1.0 



Table 1: Parameters in our Monte Carlo study. This sequence corresponds to the fixed physical 
lattice size Lg = 4.0. 



the lattice Dirac operator D by using the LU decomposition. We do not introduce any 
(super symmetry breaking) mass terms of fermions and bosons. 

We carried out simulations with the parameters in table |]. The sequence, according to 
eqs. ( |2.32 ) and ( |2.33| ), corresponds to a fixed physical lattice size Lg = 4.0 and the lattice 



spacings ag = 1.0, 0.8, 0.666, 0.571 and 0.5, respectively. For each value of (3, we stored 
1,000-10,000 independent configurations extracted from 10 6 trajectories of the molecular 
dynamics. The constant e in eq. ( p. 12 ) is kept fixed at e = 2. 6. 12 



Expressing the determinant of the Dirac operator in the form 

det{L>} = re ie , -tt < 9 < tt, (3.3) 

(generally the determinant is complex due to lattice artifacts in the present formulation) 
we evaluate the pfaffian by 

Pf{D} = ^r~e i6/2 , (3.4) 

because (Pf{-D}) 2 = det{D}. This prescription, however, may give a wrong sign for the 
pfaffian. For example, if Pf{D} = y/re^ 2 /^ m , we have = —2n/3 and the prescription (^ 
gives Pf{D} = -^/re^ 1 / 3 - 17 ™ = —^/re^ 2 ^ 3 ^ 1 which is wrong in sign. The prescription (^ 
gives the correct sign of the pfaffian, provided that — tt/2 < Arg(Pf{D}) < ir/2 (and oth- 
erwise the prescription gives a wrong sign). Although this is expected to be the case for 
large (3 (i.e., when close to the continuum), to determine the true sign of the pfaffian, we 
have to compute the pfaffian itself in some direct way. This is quite time-consuming 13 and 
we do not adopt this method in this paper. Instead, to have an idea how the prescrip- 
tion ( |3.4p works in practice, we measured the distribution of the pfaffians over a subset of 
our ensemble in table [l]. The behavior in the figure [l] clearly accords with our expectation. 
For large (3 (i.e., close to the continuum), the distribution gathers around the positive 
side of the real axis and the condition —tt/2 < Arg(Pf{L>}) < 7r/2 is fulfilled. Even for 
the smallest (3 in our simulation, (3 = 4.0, the distribution is significantly biased on the 
side of the positive real axis. Thus the systematic error introduced by the wrong-sign de- 
termination due to the prescription ( |3.4[) would be negligible compared to the statistical 
error. 



12 We observed a tendency such that the autocorrelation time becomes shorter for smaller e. Thus small e 
would be favorable from a viewpoint to accumulate a large number of configurations. On the other hand, 
it appears that smaller e implies smaller fluctuation of distribution of configurations and might be disad- 
vantageous from a viewpoint of the overlap problem. We did not systematically investigate this problem of 



an optimal e. Our present e is rather large in view of eq. (2.15). 

13 It can be seen that the algorithm for the pfaffian (appearing, for example, in ref. Q|) is an 0(n 4 )-process 
for a 2n x In matrix, while the LU decomposition has an 0(n 3 )-process algorithm. 
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Figure 1: The distribution of pfaffians in a subset of the quenched ensemble used in our simulation. 
The phase, Arg(Pf {£>}), and the modulus in logarithm, log 10 (10 16 | Pf{£)}|), are plotted in the polar 
coordinate. The number of samples is 1,000 and 100 for (3 = 4.0 and ft = 16.0, respectively. 



We consider also the quenched approximation, i.e., 

((^)) quenched = ^Tp ( 3 - 5 ) 

This provides a useful standard with which one can observe the extent of the quantum 
effect of dynamical fermions. 



3.2 One-point supersymmetric WT identities 

First, we consider supersymmetric one-point WT identities implied by the exact Q-invariance 
of the lattice action. These are of the form (((^(something))} and identically vanish because 
of Q-invariance of the action and the integration measure. These should hold for any lattice 
parameter, if the integration (especially that over fermionic fields) is properly performed. 
Thus, from their validity in numerical simulation, we can confirm the correctness of our 
code/algorithm. In particular, we can observe whether the re- weighting method works or 
not. 
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Figure 2: Expectation values of $Z i=1 £b»(:e) — (3/2) (JV^ — l)a 



Since the lattice action (|2 
density, 



is Q-exact, we have ((5)) = 0, or, in terms of the action 



tr\H(x) - -i*TL(a?) 




0. 



(3-6) 



One may further simplify this relation. The second term is the expectation value of the 
action density of the fermionic fields. Since the action is bi-linear in fermionic fields, we 
have J2t=l i( £ Fi(x)}) = -2{N 2 - l)a~ 2 (the coefficient 2 is (1/2) x 4, where 1/2 reflects 
the Majorana nature of the system and 4 is the number of fermion species). Similarly, the 
auxiliary field H(x) can be integrated out and the last term becomes (1/2) (iV 2 — l)a~ 2 
after integration. Thus 

3 

^((£ Bi (x)))-|(iV c 2 - 1)1 = 0. (3.7) 
£ a 

i=i 

In figure we plotted the left-hand side of this relation (in a unit of g 2 ) as a function of the 



lattice spacing ag. The real part is consistent with the expected identity (3.7) within la 
for all values of ag, except ag = 0.571 that is 1.5a away. This agreement strongly indicates 
the correctness of our code/algorithm. The average of the imaginary part is consistent 
with zero, as it should be, although its fluctuation is comparable to that of the real part. 

What is intriguing with figure § is that one can see clear distinction between the 
re- weighted average (|3.2|) and the quenched average ( plq) . This illustrates that the re- 
weighting method works very well and the effect of dynamical fermions is properly included 
(at least for the present quantity). A perturbative argument shows that each term of the 
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Figure 3: Expectation values of Cbi(x) + Cfi(x). 



action density behaves as 

«iCB2(x)» ~ \{N 2 C - 1)1 ((Cb3(x))) ~ (N 2 - 1)1, (3.8) 

in a — > because of one-loop diagrams and ((jCbi(^))) starts with a two-loop diagram which 
behaves as ~ 

(ln(a/L))V- The leading 0(a~ 2 ) singularities are thus cancelled out in the 
sum X^=i ((^Bi(^))) — (3/2) (iVp — l)a -2 and this leaves a function of the form f(a/L, Lg)g 2 . 
This function identically vanishes if supersymmetry holds, but it is a non-trivial function 
in the quenched approximation. What is shown in the figure ^ with "quenched" is this 
function. 14 

The identity ( |3.6[ ) can be divided into several pieces, each of which should hold sepa- 
rately. The first one is 



«Q0i(x)» = ((C B i(x))) + «£fi(*)» = 0, (3.9) 

and the left-hand side is plotted in figure ||. The relation is confirmed within 1.5a expect the 
case ag = 0.571. Note the difference in scale of vertical axis compared to figure H Although 
the results with a quenched ensemble are certainly inconsistent with the supersymmetric 
relation ( |3.9[) , we do not see clear separation between the re- weighted average and the 
quenched average. This seems to be related to the fact that diagrams that contribute to 

14 The lattice perturbation theory is not useful to evaluate this function even for ag — > 0. In the loop 
expansion, we may have terms of the form, say, (\n{a/L)) e (Lg) 2 ( e ~ 2 ^ g 2 at ^-loop level. All this type of 
terms equally contribute to the function / in the present continuum limit, in which Lg is fixed (recall that 
Lg = 4.0 in our simulation). 
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Figure 4: Expectation values of 2Cb2(x) + Cf5(x). 

((jCbi(x)}} and ((£fi(x))) and contain virtual fermion loops start with three loops, a rather 
higher order. 

Another piece of eq. ([TB) is 



«Q0 2 (*)» = -L ({tr{-iH(x)<t> TL (x)})) + ((£ F5 (x))) = 0. 



(3.10) 



Under the gaussian integration, the auxiliary field can be replaced by H(x) = ^$tl(^) 
and the above becomes 



2((£ B2 (x))) + ((£f5(x))) = 0. 



(3.11) 



In figure [|, the left-hand side of this relation is plotted. The global feature is similar to 
that of figure ^ and the relation is reproduced within almost la. 

The situation is again similar with the last piece of the relation (|3.6| ): 



((Q0 3 (x))) = ((C B3 (x))) + ((C F3 (x))) + «£ F 4(x)» + «£f 6 (x)» = 0, (3.12) 

whose left-hand side is plotted in figure || 

So far, we have observed WT identities implied by the exact Q-symmetry of the lattice 
action. The continuum theory fl2.2|), on the other hand, is invariant under also other 
fermionic transformations, Qoi, Qo and Q\. In the lattice framework, the invariance under 
these transformations is expected to be restored only in the continuum limit. The fermionic 
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Figure 5: Expectation values of Cb3(x) + Cf3(x) + £f4{x) + Cm{x) 



transformation Qoi is given by 

Qmcj) = 0, 
QoiV = 2H, 



QmH = - [(f), i]}, 



Q 



01< 



-2x, 



l, 



that can be obtained by following substitutions in the Q-transformation (J2T 

1 1 



2 V 



-x, 



(3.13) 



(3.14) 



where eoi = — eio = 1- Since the action ( [2.2| ) is invariant under these substitutions, the 



invariance of the continuum action under eq. ( 3.13 ) is obvious. Associated with this Qoi- 
invariance, in the super symmetric continuum theory, we have 



Qoi^ tr j-^x[$ 



continuum 



+ ±((tr{- x [<f>, x ]})) 



continuum 9 

Thus, corresponding to this relation, one might expect 

((£bi(x))) + ((£ F2 (x)))-0 
holds in the continuum limit a — > 0. 



continuum 



0. 



(3.15) 



(3.16) 
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Figure 6: Expectation values of Cbi(x) + Cf2{x). 



In figure ||, we plotted the left-hand side of eq. ( |3.16| ). It appears that the average 
approaches a non-zero number around 0.15, instead of zero (the imaginary part is consistent 
with zero, as it should be). This does not contradict with the supersymmetry restoration. 



As already noted, the argument [15] for a restoration of supersymmetry in the continuum 
limit is not applied to correlation functions containing composite operators. In particular, 
there is no general guarantee that the bare WT identity ( |3.16| ) holds in the continuum 
limit. 

We note that if supersymmetry in the 1PI effective action is restored in the continuum 
limit, it is UV finite, that is, all 1PI diagrams with external lines of elementary fields 
are UV finite. Power counting (taking gauge invariance into account) shows that only 
scalar mass terms suffer from superficial UV divergence. Scalar mass terms are, however, 
inconsistent with supersymmetry. So if the 1PI effective action is supersymmetric, it is UV 
finite. 15 On the other hand, composite operators Cbi(x) and Cp2(x) induce logarithmic UV 
divergence at two-loop level. If supersymmetry of the 1PI effective action is restored, this 
two-loop level divergence, caused by the presence of composite operators, is the only source 
of UV divergence in {(Cbi(x))) and ((£f2 (#)))• Moreover, that remaining two-loop level 
divergence is cancelled out in the sum ((Cbi(x)}} + ((Cf2(x)))- This argument shows that, if 
supersymmetry in the 1PI effective action restores, the dependence of ((£bi(x)}} + ((Cf2(x)}} 
on ag decreases as ag — > 0, i.e., it approaches a constant (but not necessarily zero). The 
behavior in figure ^ is consistent with this picture based on a restoration of supersymmetry. 

What is not completely clear to us is that even the quenched average seems to have 



The contrary is not true. The UV finiteness of the effective action does not imply supersymmetry, as 
finite scalar mass terms are allowed for the former. 
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the same behavior. Actually, within almost la the re-weighted average and the quenched 
average are degenerate. So, although figure |6| is consistent with a scenario of a supersym- 
metry restoration, we cannot conclude the restoration of supersymmetry from the above 
result. 

The continuum action is invariant under also 



QqAq 

Qo(t> = 
Qo4> = 



l 



-X, 



-H, 



Qov = 
Qox = 

QOH: 

Qoipo 



-2iD </> : 
iDi<f>, 

~- [0,^1], 

= -0,4 



that can be obtained by the substitutions in eq. (2.3) 
1 , 1 



2 V 



^0, 



X 



2* 



-x, 



(3.17) 



(3.18) 



Corresponding to this symmetry, we have 



Qo-ntr 



r 



continuum 



1 



-2 tr \-A^\ 



9 



and one might expect 



continuum 



+ -2 <M-^hta0}» 



continuum 



(3.19) 
(3.20) 



((Cm(x))) + «£f3(z)» -> 0, 

in the continuum limit a — > 0. The result (figure 0) is similar to the previous one. The 
average seems to approach a non-zero number around 0.05 and we may repeat the above 
argument. 

Another fermionic symmetry Q\ is obtained by further exchange tpo ip\ in eq. ( 3.1S| ). 
Corresponding to this, one might expect 



((£bi(x))) + ((£ F 4(x)))^0. 



(3.21) 



The result of numerical study is plotted in figure [8|. The result is very similar to that of 
figure 0. 

3.3 Expectation value of scalar bi-linear operators 

To illustrate possible use of lattice simulation of the present kind, in this section, we consider 
expectation values of gauge-invariant bi-linear operators of scalar fields, a~ 2 tr{(p(x)<p(x)} 
and a~ 2 ti{<f>(x)(f)(x)} (the factor a~ 2 is multiplied for the rescaling ( 2.34| )). The classical 
action of the 2d N = (2,2) SYM vanishes identically for all configurations of constant 
scalar fields such that [4>,(j>} = (other fields are set to zero). These are so-called flat- 
directions and classical vacua are infinitely degenerate. Moreover, this degeneracy is not 
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Figure 7: Expectation values of Cbi{x) + £f3(x). 
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Figure 8: Expectation values of £bi{%) + £f4,(%)- 



lifted by radiative corrections to all order of perturbative theory. Thus, the expectation 
values of scalar fields in quantum theory are of interest and, if Monte Carlo simulation is 
useful, a prediction on these expectation values should be feasible. 

First, we consider a~ 2 tr{(f>(x)(f)(x)}. This operator is invariant under the global U(l)^ 
transformation, which acts on scalar fields as (f)(x) — ► e 2ia (p(x) and 4>{x) — > e~ 2ia 4>{x) . 



-17- 



L/a = N 


8 


7 


6 


5 


4 


c(a/L) 


0.379295 


0.357928 


0.333234 


0.304 


0.268229 



Table 2: The counter constant c(a/L) given by eq. (3.23) 



The continuum limit of this quantity itself is meaningless, because it is a bare quantity 
and suffers from UV divergence. It should be renormalized. A power counting argument 
shows that the superficial UV divergence comes from the simplest one-loop diagram and 
the divergence is logarithmic ~ ln(a/L)g 2 . If supersymmetry of the 1PI effective action 
is restored in the continuum limit, as we assume at the moment, this one- loop divergence 
is the only source of UV divergence of ((a~ 2 ti{4>(x)(j){x)})) (recall the argument below 
eq. (p6D ). 

So we define the renormalized operator (the normal product) 

M[a- 2 tr{0(x)0(x)}] = a' 2 tr{(j)(x)0(x)} - (N 2 - l)c{a/L)g 2 (3.22) 

by subtracting a c-number, the value of the one-loop diagram. This subtraction must 
remove all the UV divergence of the composite operator. This simplicity is a special 
property of the present two-dimensional (supersymmetric) theory. 

The coefficient c(a/L) of the counter constant is given by a simple scalar one-loop 
diagram and, on a finite size lattice, it is 

N-l N-l 

c(a/L = l/N) = — 2 £ £ _ — — . (3.23) 



rao=0 rti=0 



At=0 



2tt 
cos — nfl 



As possible prescription for the zero mode, we do not include (no,m) = (0,0) in the sum. 
Values of this counter constant are listed in table || for the cases in our simulation. 

The result of our Monte Carlo simulation is figure ||. 16 First of all, we see clear 
separation between the re- weighted average and the quenched one. The difference is thus 
due to the effect of dynamical fermions. This effect uplifts the expectation value and this is 
consistent with the picture that, in quenched (i.e., non-supersymmetric) theory, the scalar 
potential is lifted by radiative corrections, suppressing quantum fluctuation of scalar fields. 
As discussed for the WT identity ( 3.16j ), if the supersymmetry is restored in the continuum 



limit, the expectation value (^J\f[a~ 2 tr{0(x)0(x)}])) is expected to become independent of 
ag as a — > 0. The behavior in figure || is more or less consistent with this expectation, 
although clearly we need further data at smaller values of ag to conclude this. In any case, 
interestingly, the expectation value appears to approach some finite number (in a unit of g 2 ) 
in the continuum limit after the renormalization ( |3.22 ). (Without the renormalization (the 



subtraction), there is a tendency that the expectation values grow as a — > 0.) The limiting 
value of ((N[a~ 2 tr{</>(a;)</>(x)}])) at a — > (while fixing Lg) in the figure itself has no 
direct physical meaning because it can freely be shifted by a further finite renormalization. 



16 We confirmed that the imaginary part is almost negligible (as it should be) and it is not plotted in the 
figure. 
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Figure 9: Expectation values of Af[a 2 tr{(j>(x)(j>(x)}}. 




-0.6 



Figure 10: Expectation values of a 2 tr{4>(x)(j)(x)}. 



However, the limiting value should depend on Lg and this dependence can be a non-trivial 
prediction. We need a much finer lattice, of course, for an extrapolation to the continuum. 

In figure 10, we have plotted (^a^ 2 tx{(p(x)(j)(x)})). For this, a perturbative argu- 
ment indicates that there is no need of renormalization. The result is clearly shows 
((a -2 tic {(ft (x) (ft (x)}^ ~ 0. This might be suggested from the fact that in two dimen- 
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sions the global U(l)# symmetry cannot be spontaneously broken, although this argument 
is not rigorous because we are studying a system in finite volume. 



4. Conclusion 

In this paper, we presented the results of our preliminary numerical study of Sugino's lattice 
formulation of the 2d M = (2, 2) SYM. By confirming WT identities associated with an 
exact fermionic symmetry of the formulation in fair accuracy, we infer that the re- weighting 
method for the dynamical fermions basically works in this two-dimensional system. On the 
other hand, although we could not conclude the restoration of full supersymmetry from 
the numerical results, all the results are consistent with the basic idea of a supersymmetry 
restoration. We computed also the expectation values of scalar bi-linear operators to 
illustrate the usefulness of this kind of lattice simulation. 

In this paper, we did not try to measure any two-point correlation function or ex- 
tended observables like Wilson loops, because it is clear that our lattice is too small to 
extract any useful information from such quantities. Interesting physics of this system is, 
of course, contained in these observables. For example, the most direct way to examine 
the restoration (and/or the spontaneous breaking) of supersymmetry is to study the mass 
spectra and two-point functions containing the supersymmetric current. For an interesting 
property of a two-point function that contains the U(1)r current, see ref. Thanks 
to FermiQCD/MDP |33|], our code is executable also on a large PC cluster without any 
change. Having obtained encouraging results in this paper, in the near future, we hope to 
report results of full-scale simulation using much larger lattice. 

There exists a natural generalization of the present manifestly Q-invariant lattice for- 
mulation to the 2d N = (4, 4) SYM [@, § and to the 2d M = (8, 8) SYM (the second paper 
of ref. [p^l ). The latter theory is especially of interest as an effective theory that describes 
the dynamics of Dl-brane. We do not find any real difficulty to set up the corresponding 
Monte Carlo simulation similar to that of the present paper. This is an interesting future 
problem. 
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